Wavelets as basis functions in canonical quantization 
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Canonical quantization of electromagnetic field is traditionally done using plane waves, ft is 
possible to formulate the quantization using other complete set of basis functions. Wavelets are a 
special kind of functions which are localized in real as well as in Fourier space. In this paper we show 
how wavelets can be used as basis functions in canonical quantization. A countable set of mode 
functions are obtained. The general formalism of the change of basis is the same for all wavelets 
which satisfy a multiresolution analysis. 
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I. INTRODUCTION 



Canonical quantization of electromagnetic field is tra- 
ditionally done using plane waves. The field is enclosed 
into a cubic cavity and the field operators are expanded 
using eigenfunctions of the cavity, i.e., plane waves. After 
quantization it is possible to take the limit of infinite cav- 
ity and eliminate the unphysical finite size cavity. The 
inconvenience of the use of plane waves is that they are 
delocalized in real space. This makes it difficult to for- 
mulate for example photodetection theories. Photode- 
tectors measure the field locally and are sensitive over a 
finite bandwidth of frequencies. It is possible to use any 
complete set of basis functions in the quantization , 
so a formulation with a more localized basis functions 
would be desirable. 

The theory of multiresolution analysis (MRA) has 
been under intensive study during the recent years ■ 
A complete set of basis functions in MRA are called 
wavelets. Wavelets are localized in real as well as in 
Fourier space. They are parameterized by scale and 
translation parameters which both get integer values. 
The translation parameter translates the wavelet in real 
space. For large negative values of the scale parameter 
the wavelet is wide and for large scale parameter values 
narrow. Typical wavelets are orthogonal with respect 
to both indices. There are many wavelets with differ- 
ent characteristics. Differences include how well wavelets 
are localized, what kind of Fourier transforms they have, 
whether they are real or complex and how symmetric 
they are. It is interesting that some wavelets have com- 
pact support, i.e., they are zero outside a certain finite 
length interval. These kind of wavelets do not have an- 
alytical expressions. The main theory of multiresolution 
analysis is the same for all useful wavelets. 

In this paper we show how wavelets can be used as 
basis functions in the canonical quantization. Real and 
orthonormal wavelets are used. New mode functions and 
operators are linear transforms of plane waves and the 
corresponding operators. Different vector valued mode 
functions for electric and magnetic fields are obtained. 
New creation and annihilation operators are the same 



for both fields and satisfy bosonic commutation relations. 
This means that formalism of the new operators remains 
the same and makes it easy to use the new basis. New 
mode functions are localized and have similar properties 
as wavelets, for examble they are parameterized by scale 
and translation parameters. 

In Sec. II we give a short introduction to the theory of 
multiresolution analysis, scaling functions and wavelets. 
In Sec. Ill we derive equations for field operators in 
wavelet basis. In Sec. IV the theory developed in earlier 
sections is applied to some simple quantum mechanical 
simulations and the wavelet and plane wave bases are 
compared. Finally in Sec. V we give our conclusions and 
suggest several generalizations of the theory developed in 
this paper. 



II. BASIC PROPERTIES OF SCALING 
FUNCTIONS AND WAVELETS 

A. Multiresolution analysis and wavelets 

In the following we give a brief introduction to the 
basic properties of wavelets. The discussion follows 
books 1^-^. We start with the scaling function 4>i{t) = 
— I & Z, (j) G L^, which spans the function space 

Ao 

/(t) = ^az0z(t), f{t)eAo. (2.1) 
I 

Introducing a new index s by the formula 

= 2^/''/'(2'i-0, s,leZ (2.2) 

we can define different function spaces Ag. If s is large, 
s 3> 0, the scaling function is narrow and peaked around 
some center point. Large negative values s <C gives a 
function which is wide. The parameter k translates the 
scaling function in a given scale by a factor 2^". 

In order to the scaling function to satisfy a multires- 
olution analysis (MRA) the function space Ag+i must 
include the function space with a lower index As 



As c A, 



S+17 



se z 



(2.3) 
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When the parameter s approaches infinity we get a 
space of square integrable functions A^o = and when 
s approaches minus infinity the result is a null space 
A-ao = {0}. Because of multiresolution analysis, it is 
possible to expand every function in a function space A^ 
using basis functions of function space Ag+i. Specifically 
we can expand scaling functions (t>s,i{t) in As using basis 
functions of As+i- This gives 



oo oo 



(j){t) = ^h{n)V2(j)(2t-n), n€Z 



(2.4) 



with some coe fficients h{n). We have chosen s = and 
used Eq. 

The wavelet space is defined to be a difference space 
between different function spaces spanned by the scaling 
functions. We define IF, as 



A 



s+l 



As®Ws 



(2.5) 



Intuitively Wg contains functions which must be added 
to As in order to get Ag+i- The basis functions in the 
function space Wg are called wavelets ■i/'(0- functions 
in Ws can be expanded using translations of a fundamen- 
tal wavelet, which are obtained using the same equation 
as for scaling functions, Eq. (|2.2|). Because Wg C A^+i, 
the wavelet function in Wg can be expanded using scaling 
functions in Ag+i. We get 



hi{n)V2(t){2t-n), n€Z, 



(2.6) 



which corresponds to the equation ( 2A ) for i/i. The coeffi- 
cients h{n) and hi{n) are not independent. They satisfy 
the relation 



hi{n) = (-l)"/i(l-n). 



(2.7) 



Using Eq. (2^) in Eq. (2_^) wc can decompose the 
right hand side into subspaces with lower indices. After 
one iteration we get 



A. 



+2 



A. 



+1 



+1 



Ag®Wg® Wg- 



(2.8) 



We can do the iteration repeatedly and in the infinite 
limit we get 



^Ag®Wg® W, 



s+l 



(2.9) 



where A^o — has been used. This means that every 
function which belongs to can be expanded as 



f{t)= J2 CgA2'^^<Pi2H-l) (2.10) 

OO OO 

J2 E dg,,ir'^^^{r't-i). 



1 — — 00 



s'—s 1— — 00 



The parameter s is totally arbitrary. We can take a limit 
s — > —oo and because A-^o = {0} the scaling function 
space can be eliminated giving 



s — — oo — oo 
oo oo 

= E E dg,ir'^i>{rt~ii (2.11) 



s — — oo — oo 



i.e., any L^-function can be expanded using translated 
and scaled wavelets. In this paper we use only real or- 
thonormal wavelets which have the property 



'ijjgi{t)ipg>i'{t)dt = Sgg'Sii'. 



(2.12) 



Multiplying Eq. ( |2.11 ) by parts with ■ipgiii{t) and inte- 

oo 

grating J dt wc get the coefficients 



(2.13) 



dg,i = / f{t)tl^g,i{t)dt. 



It follows from the multiresolution analysis that the in- 
tegral over the wavelet gives zero, i.e. 



'4}{t)dt = 0. 



(2.14) 



This means that wavelets must have some kind of oscil- 
lating structure as the name suggests. 

So far we have not defined any wavelets explicitly. 
There are several methods to construct wavelets for dif- 
ferent purposes. One method is to divide the Fourier 
space in such a way that the resulting wavelets fulfill 
multiresolution analysis requirements and are orthogo- 
nal. The wavelets obtained using this method can be well 
localized but do not have a compact support. Another 
method is to derive wavelets based on the filter coeffi- 
cients. The filter coefficients of typical wavelets satisfy 
the fundamental condition 



1 

7i 



and are orthogonal 



Y,h{n)h{n-2k) = 5k^. 



(2.15) 



(2.16) 



If the length of the filter coefficient sequence is long 
enough these two conditions do not determine the co- 
efficients uniquely. The remaining degrees of freedom 
can be used to give additional desirable properties to 
wavelets. Typical choices are to demand the wavelet or 
scaling function to be as smooth as possible. These kind 
of wavelets can have a compact support, i.e., they are 
zero outside a specific region. 
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B. A few examples of wavelets 

One of the best known scaling function is a sine- 
function 



(x) = sinc(a;) 



sin(7rx) 



(2.17) 



with sinc(O) = 1. The corresponding wavelet is ipi^) — 
2(f){2x) — (jii^)- Figure |^ shows the scaling function and 
the wavelet which are both localized around a; = 0. Both 
functions are not so well localized as would be desirable 
for many practical applications. The Fourier transforms 
of the scaling and wavelet functions are 



0H = 



0, otherwise 
0, otherwise. 



(2.18) 
(2.19) 



The reason why oscillations do not decay rapidly is be- 
cause of the abrupt changes of the Fourier transforms 
at |w| = TT and |a;| = 27r. One specialty of the sinc- 
family of wavelets is that the division of frequency space 
to different scales is orthogonal, i.e., there is no overlap 
of the Fourier-transforms of wavelets with different scale 
parameters s. Because of poor localization in real space, 
Shannon wavelets are rarely used. 



Scaling function 



Wavelet function 
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FIG. 1. Fundamental Shannon scaling and wavelet func- 
tions {(j) and ■(/)). Both functions have analytical expressions, 
are symmetric and have an oscillating structure. Oscillations 
are large also far away from the origin. 

It is possible to smoothen the change in Fourier trans- 
forms in such a way that orthogonality and multiresolu- 
tion analysis requirements are preserved. The resulting 
wavelets are called Meyer-wavelets. There are several 
different Meyer scaling and wavelet functions depending 
on the smoothening function. The Meyer scaling and 
wavelet functions used in this paper are shown in Fig. ||. 
They are clearly more localized than the corresponding 
Shannon functions. Absolute values of the Fourier trans- 
forms of scaling and wavelet functions are shown in Fig. 
^. The Fourier transform of a scaling function (f> is flat 
around lo = 0. Around frequencies uj ~ ±tt the Fourier 



transform decays to zero and is zero at larger values. 
The smoothing compared to Shannon case is clearly seen 
in both functions. For wavelets two different scales are 
shown. As parameter s increases the transform is wider 
and shifted further away from the origin. 
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FIG. 2. Fundamental scaling and wavelet functions of 
Meyer wavelets. Both functions are better localized than the 
corresponding Shannon functions (Fig. Unlike the Shan- 
non functions, the Meyer scaling function and wavelet are not 
totally symmetric. 
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FIG. 3. Left: The absolute value of the Fourier transform 
of a Meyer scaling function. The transform is nonzero at small 
frequencies and decays to zero at around |fe| = n. The change 
is smoothened in order to improve localization. Right: The 
absolute value of the Fourier transform of a Meyer wavelet at 
two different scale parameter values s = and s = 1. There 
are two peaks symmetrically around zero. As the scale pa- 
rameter increases the peaks become lower and are shifted to 
higher frequencies. 



Th e len gth of filter sequences h{n) and hi{n) in Eq. 
( ^.4|) and (2^) for sine and Meyer scaling fuctions and 
wavelets are infinite. As a result of this wavelets do not 
have a compact support but are nonzero, although pos- 
sibly very small, over the whole x-axis. It is possible to 
construct wavelets which have compact support, i.e., they 
are exactly zero, not just exponentially small, outside a 
specific region. The length of the filter sequence for these 
kind of wavelets is finite and the corresponding wavelets 
and scaling functions do not have analytical forms. One 
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group of these wavelets are the Daubechies wavelets. In 
this group of wavelets there are several wavelet families 
which can be characterized by the lengths of filter se- 
quences. As the length of the sequence increases the 
wavelets and scaling functions become smoother. Figure 
^ shows Daubechies scaling and wavelet functions with 
a filter coefficient length = 6 in Eq. (2.2) and (2.6). 
Both functions are nonzero in the interval < a; < 5. 
They are continuous and derivable. The functions are 
not symmetric. The general behavior is that the scaling 
function is concentrated to the left of the nonzero inter- 
val and the wavelet is significantly nonzero in the middle. 



Scaling function 



Wavelet function 





FIG. 4. Daubechies scaling function and wavelet with a 
filter length A*' = 6. Both functions are continuous and deriv- 
able. They have a compact support, i.e., they are strictly zero 
outside an interval < a; < 5. 



C. Wavelets in higher dimensions 



the same in both directions. Note that the scaling of 
the wavelet is different compared to the one dimensional 
case. Any two-dimensional L^-function can be expanded 

as 



/(r)=^4V4i(r)^ E E E E^^i^^iW- 



sli 



-co l^— — oo ly — — oo i—1 



(2.23) 



Shannon and Meyer wavelets were constructed by di- 
viding the Fourier space. For these wavelets the three dif- 
ferent wavelets in two dimensions, obtained by multiply- 
ing one dimensional scaling functions and wavelets, are 
significantly nonzero in different regions of k-space. Fig- 
ure 1^ shows nonzero regions for three different wavelets. 
Wavelets ip^{r) and ip'^{r) are nonzero in regions \kx\ < 

^mirij ^min — \^y\ — ^max and \ky\ ^ ^mm ; ^min — 

\kx\ < kmax respectively. Wavelet V''^(r) is nonzero in 



"min ^ I ky I 



^ kri 



With 



regions kmin ^ \kx\ ^ kmax-, k 

a smaller scaling index s, the scaling function in the cen- 
ter is split into similar regions with different wavelets. 
With larger parameter s, larger k-values are divided to 
regions. The division of frequency space is exact only for 
sine- wave lets. For all other wavelets the Fourier trans- 
form has nonzero values also outside of its main region. 
How well the Fourier transform is localized into the re- 
gions depends on the type of a wavelet used. For wavelets 
with compact support and small filter sequence length 
the division described above is not very well. 



The discussion above has been about one dimensional 
wavelets. We want to use wavelets as basis functions 
in canonical quantization so three dimensional wavelets 
must be used. One method to construct multidimen- 
sional wavelets is to use products of one dimensional scal- 
ing and wavelet functions 1^,0]. The following products 
are all wavelets in two dimensions 

i}\r) = 0(a;)V'(2/) 

V''(r)=^(x)0(y) (2.20) 
V'^(r) = il;{x)ipiy). 

The integral over all functions above vanishes as is re- 
quired for wavelets. The two dimensional scaling function 
is a product of two one dimensional scaling functions 

0(r) = (f>{x)(f>{y). (2.21) 

The wavelets with different scaling and translation pa- 
rameters become 

^:,(r) = Vjo(r - 2-^1) = 2«^^o(2^r " !)■ (2-22) 

The translation parameter 1 is now a vector with two 
components which give the translation in x and y direc- 
tions. The scaling parameter is a scalar, so scaling is 



k 



¥(k,)V(ky) 




V(kjvii(ky) 




<l'(k,)(Kkj) 


V(kJ*(kj) 


V(k^)>|/(kj) 


*(kjv|/(ky) 


VCkJvCkj) 


-k 
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FIG. 5. Regions of fc-space where a two dimensional scal- 
ing function and different wavelets are nonzero at a specific 
scale. The scaling function is nonzero at small frequencies. 
Fourier transforms of ip^ = 4>{kx)ip{ky) and ip'^ = il){kx)4>{ky) 
are nonzero at two different regions around ky and k^ axes. 
The Fourier transform of the third wavelet ^/)^ = ip{kx)4>iky) 
is nonzero at the corners of the frequency interval. 
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Three dimensional wavelets can be constructed in a 
similar way as two dimensional ones. We get one scaling 
function and seven wavelet functions which are products 
of one dimensional scaling and wavelet functions 



0(r) = (j){x)(l){y)(j){z) 
^2(r) = 0(a:)V(y)0(^) 
^4(r)=V(x)</.(j/)0(^) 
V'^(r)=V(x)V'(y)0(^) 



V'^(r) = V(a;)(/>(j/)V'(^) 



(2.24) 



The division of k-space is a direct generalization of the 
two dimensional case. The scaling relation becomes 

V'li(r) = rsoir - S-'^l) - 23^/2V;S„(2^r - 1). (2.25) 



The expansion of L^-function is given by Eq. (|2.23| ) with 
the difference that the 1- vector is now three dimensional 
and index i gets values from one to seven. 

D. Expansion of a plane wave using wavelets 

In the following section we need expansions of plane 
waves using wavelets. The normalized plane wave in a 
box of side L can be expanded as 



L3/ 



^e^''-- = ^dL>:,(r). 



(2.26) 



sli 



Multiplying by parts with (r) and using orthogonal- 
ity of wavelets gives the coefficients 



"k,sl — ^3/2 



(2.27) 



The factorization is possible because multidimensional 
wavelets are products of one dimensional scaling func- 
tions and wavelets. The index ix,y,z is either or 1 de- 
noting scaling and wavelet functio ns re spectively. Simi- 
l ar fa ctorization is possible in Eq. ( ^.26 ). Using equation 
( ^.25 ) we get after the change of integration variable 



- exp(z2-^k . l)dl,, = 2-3^/2 exp(z2-^k • l)4-.k,oo- 

(2.28) 

Thus every coefficient dj^. can be calculated using the 
translated and scaled Fourier transform of a fundamen- 
tal wavelet. The scaling parameter s changes also the 
absolute value of the coefficients whereas the translation 
parameter 1 gives only a phase shift. Because of orthog- 
onality of plane waves and wavelets we get the following 
relations 



k,sl"k,s'r 



Sss'Sll'Sr, 



(2.29) 
(2.30) 



Next we study the parameters dj^. more closely. We 
restrict the discussion to one dimension, i.e., the param- 
eters k and 1 are one dimensional. Figure ^ shows the 
absolute values of the coefficients with several different 
fc-valu es for Meyer wavelets. Note that according to Eq. 
( ^.27 ) the translation parameter gives only a phase shift. 
The coefficients of a plane wave are shown with three 
different fc-values. In all cases coefficients with only a 
few scale parameters are nonzero. For small fc- values 
k = 0.39, two scales with parameters s < have nonzero 
values. Remember that the Fourier transform of the fun- 
damental (s — 0) Meyer wavelet is mainly nonzero in the 
interval tt < k < 27r. The fc-value k = 4.71 = I.Stt is care- 
fully chosen to be in the middle of the scale s — inter- 
val. It is seen that coefficients for all other s-parameters 
are zero. With large fc-values the peaks are centered at a 
higher scale as can be seen from the last peak. Peaks with 
large scale parameter are smaller than peaks at small 
scale. This is because of the factor 2^^'^^^ in equation 
( p. 28 ). One could think that at large scale the peaks are 



lower because there are more Z-values which share the 
contribution from the single plane wave. Remember that 
at high scale values the shift of a wavelet is smaller when 
the translation parameter is changed by one unit. 



k=0.39 
k=4.71 
k=100.0 




2 
scale parameter 



FIG. 6. The absolute value of the Fourier transform of 
a wavelet dk,si with three different fe-values. For k = 0.39 
scales s = —4 and s = —3 are nonzero. Because the value 
k = 4.71 = I.Stt is in the middle of the scale interval, only 
one s-value is nonzero. For k = 100.0 the two scales are 
nonzero. 
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III. CANONICAL QUANTIZATION USING 
WAVELETS 

Traditionally canonical quantization of field is done us- 
ing plane waves as basis functions ||]. The field is en- 
closed to a cubic of side L. All functions can be expanded 
using the eigenfunctions of the quantization volume, i.e., 
plane waves. The quantization itself is done by intro- 
ducing bosonic operators for every basis function. The 
vector potential of field can be expanded 



A = 



1 



i3/2 

ko- 



2wkeo 



1/2 



(Oko-ekcre 



h.c. 



(3.1) 



The vector eto- gives the polarization of the plane wave 
with wave vector k and polarization a. Operators ako- 
and a\^^ are bosonic annihilation and creation operators. 
The parameter k gets values 



2imi 27rn2 27m3 



For electric and magnetic fields we get 

1 



ni,2,3 = 0,±1,±2... (3.2) 



E(r) 



L3/2 



E 



1/2 




1/2 



E(r) 



aA(r) 



1/2 



L3/2 V 2eo ; 

kcr ^ ^ 



B(r) = V X A(r) 
1 



■E 



kCT 



1/2 



(iakcr(k X ekcr)e 



ik-r 



(3.3) 

(3.4) 
h.c.) . 



Field operators above are in plane wave basis, which is 
parameterized by k-vectors. 

Next we change the basis from plane waves to wavelets. 
We could just expand the plane waves usin g wavelets and 
use this expansion in Eqs. (3.3) and (3.4). However this 
approach would not allow the separation of mode func- 
tions and operators. Therefore we proce ed b y inserting 
a unity in the form ^ (5kk' to expansions ( |3.3|) and ([3.4|). 

k 

The sum is now divided in such a way that all other k- 
values exce pt the operator k are changed to k'. After the 
use of Eq. ( ^.29 ) for the delta function we get 



h.c. 



^ik' -r 



"2^ 



1/2 



ekvrf- 



V«k'.sl' 



iU. -r 




(3.5) 

(3.6) 

(3.7) 
(3.8) 



where the new basis functions and operators are 



L3/2 Z^\2eo J 



'^ktslfkae 



ik-r 



(3.9) 



= u-(r - 2-^) = 2-u-,(2^r - 1) (3.10) 
^k. = E'^k>ka- (3.11) 

k 

Similarly for the magnetic field we get 

= E E W + h.c.) , (3.12) 

sli (7 



where the operator ^ is given by equation (3.11) and 



<t(r) 



E 



1/2 



4%i(k X ek.)e^'^ "- (3.13) 



u:«(r-2-q) = 22Xo,.(2^r~l). 



(3.14) 



The new mode functions behave in the same way as 
wavelets when the indices are changed. The translation 
parameter 1 translates the mode functions in three di- 
mensions and the parameter s compresses and stretches 
them. Even though the quantization in the wavelet basis 
is in free space the operators and mode functions are ex- 
panded using a countable set of basis functions. On the 
contrary to plane waves the polarization of the new mode 
functions is not constant. As is the case with wavelets, 
the integral of the mode functions over quantization vol- 
ume vanishes 



j n%{v)d^v = J u:{^,(r)d3r = 0. 



(3.15) 



This means that also the mode functions have similar 
'waviness' as wavelets. 

It is easy to show that the new operators 6* j ^ and S*], ^ 
obey bosonic commutation relations 



0, 



(3.16) 
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(3.17) 



Because of this the operators act hke annihilation and 
creation operators for wavelet modes. This means that 
the formalism in wavelet basis is the same as wi th p lane 
wave operators. Multiplying both sides of Eq. ( 3.11 ) by 
'^k si performing the sum ^ we get with the help of 

sli 

the orthogonality integral 



(3.18) 



Similarly it is possible to expand plane wave basis func- 
tions using wavelet mode functions. 

The field Hamiltonian in wavelet basis is obtained as 
an integral of the energy density over the quantization 
volume 



Ht.= \ \ \t^E\v) + ^S'(r) ) rf'r, (3.19) 



where expansions dJI) and (|t|) are used for electric 
and magnetic field operators. In the calculation of ( 3.19| ) 
the following integrals are needed (last one is listed only 
for completeness) 



ZCq 

zeoc 



(3.20) 
(3.21) 
(3.22) 



where the real coupling constants of the modes are 



(3.23) 



(27r) 



(3.24) 



Using ( ^.28D we get 

<M'i' = E^k4%o4,.'oexp(ik • (2-«'l' - 2-n)) 



where 



(3.25) 



(3.26) 



transforms have a nonzero overlap. Because the wavelets 
are localized both in real and Fourier space this means 
that for majority of mode pairs the coupling is zero. The 
detailed structure of the coupling constants w^,^^, ], de - 
pends on the wavelet used. The second term in Eq. ( |3.28| ) 
is the result of the self-coupling of the modes and gives 
the zero point energy of the field. As is the case for plane 
waves the energy per mode is half of the unit excitation 
field. 



The identity ( |3.25| ) shows that the coupling constant is 
a function of the difference of t he sc aled translation pa- 
rameters only. Using equations ( ^.20 ) and ( 3.21 ) and the 
identity 

I u:f,(r) . utl,,ir)d'r = -c^ j n%„{v) ■ u;',f (r)d3r 

(3.27) 

we get after a straightforward calculation for the free field 
Hamiltonian 



(7 sli s'V- 



2 / J / J si, si" 

(7 sit 



(3.28) 



The modes are coupled with the coupling constant 

^si*s'r- '^^^ diagonal elements "w^l g\ give the energy of 
the field state which has only one mode with parameters 
s, 1 and i excite d, i.e., energy of one wavelet quantum. 
From Eq. ( |3^ ) and ( |3^ ) it is seen that the coupling 
is nonzero only if the corresponding wavelets and Fourier 



IV. NUMERICAL SIMULATIONS USING 
WAVELET MODES 

A. Conventions and wavelets used 

In this section we give examples how to apply the the- 
ory developed in the last section. We do the calculations 
in two rather than in three dimensions because visualiza- 
tion is easier and all essential features are included. Thus 
the k- vector is restricted to xy plane k = k^G^x -\- kye,i 
and the polarization vector is eki = —63. Here we 
exclude the other polarization ek2. The cros s product 
in the expansion of the magnetic field (3.4) becomes 



k X Eki — —ky&x -I- The choice of the field con- 

figuration considered is the same as in our earlier paper 
The mode functions for electric and magnetic fields 
become 



<uW = -fE(i^) <..e3e'"^ (4-1) 



<u(r) = 



(4.2) 
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where the summation over k is now over the xy plane. In 
the following we choose units such that h = eo = c = 1. 
In all examples in this section we use the Meyer wavelets 
shown in Fig. ^. Meyer wavelets were designed using 
smoothened Fourier transform of Shannon wavelets, so 
the division of the Fouricr-planc shown in Fig. ||is rather 
valid. The Fourier transforms of the mode functions are 
quite similar to the Fourier transforms of the correspond- 
ing wavelets. In all simulations the parameters are chosen 
in such a way that it is necessary to consider only two 
scales s = and s = 1. This means that the frequencies 
\kx\ < TT, \ky\ < TT and \kx\ > 47r, \ky\ > An are excluded. 
One has to note that in two dimensions the scaling of the 
wavelet mode functions become 



> 
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FIG. 8. Absolute value of magnetic wavelet mode func- 
tions u*f . The parameters are the same as in Fig. ^. Mode 
functions are localized and different from the corresponding 
electric mode functions (Fig. ^). 

B. A Gaussian photon 



U 



r 



iE.B , 



1), (4.3) 



instead of t he sc aling in three dimensions given by Eqs. 
(p^ and (pH). 



The absolute values of the electric and magnetic mode 
functions with indices s = and = ly = are shown 
in Figs. 1^ and ^. On the left in both figures there is the 
mode function with index i = 1 and on the right i = 3. 
The i — 2 mode function is the same as z = 1 case with 
axes X and y changed. All mode functions are clearly 
localized at the origin and have wavelet type of oscillat- 
ing structure. Because the Fourier transform of i = 3 
two dimensional wavelet is nonzero at higher frequencies 
than Fourier transforms of i = 1 and i = 2 wavelets, the 
oscillations of i = 3 mode function have smaller details. 
As was explained earlier the parameters s and 1 scale and 
translate mode functions. 






• • * 
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X X 

FIG. 7. Absolute value of electric wavelet mode functions 
u*f with parameters s = and 1 = 0. The function on the 
left has index i = 1 and on the right i = 3. The i = 2 mode 
function is similar to i = 1 case with x- and y-axes changed. 
Mode functions are localized at the origin and have wavelet 
type oscillations. 



In this section we study a one photon field which is 
a superposition of different single excitation plane wave 
modes. The distribution of the absolute values of the 
mode coefficients is a Gaussian centered at some point 
ko with a variance in k-space. Thus the field state is 



I*) 



k 



Ck|lk,{0})=5]ckal,|{0}) 



where 



Ck 



(27rA^ 



-1/2, 



exp 



(k-ko 
4A2 



(4.4) 



(4.5) 



The transformation to wavelet basis can be done using 
the expansion of operators in a wavelet basis (3.11). 
The state vector in wavelet basis becomes 



sli 



where 



k 



k.sl'^k- 



(4.6) 



(4.7) 



Because coefficients d^. and Ck factorize, it is possible 



to factorize the sum (4.7) 



First we briefly study a one dimensional Gaussian dis- 
tribution. We have a Gaussian photon with parameters 
ko = 18.0, xq = 10.0 and A| = 10.0. The width of the 
distribution is relatively large so one would expect several 
scales to have nonzero coefficients. Figure |^ shows coeffi- 
cients of three different scales as a function of a transla- 
tion parameter. The peak on the left has the scale param- 
eter s = 1. When the translation parameter is changed 
by one, the wavelet and mode functions are translated by 
2~'^ in real space. At scale s = 1 the translation unit is 
2~i — 0.5. The distribution is centered at xq = 10.0 so 
coefficients at s = 1 are centered at I = 20. When s = 2 
and s = 3 the coefficients around I = 40 and I = 80 have 
nonzero values. It is seen that a Gaussian state with the 
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parameters used has the biggest contribution at a scale 
s = 2. 



s=1 
s=2 
s=3 



40 60 
translation parameter 



FIG. 9. Absolute values of wavelet mode coefficients for a 
one dimensional Gaussian photon as a function of a transla- 
tion parameter. The parameters are ro — 10.0, ko = 18.0 and 
— 10.0. Coefficients are zero only for three scales. The 
peaks are at different positions because the translation unit is 
not the same at different scales. Scale s — 2 gives the biggest 
contribution. 

Next we study the time evolution of a two dimensional 
Gaussian photon. Figure |l^ shows the time evolution of 
the energy density distribution of a phot on at two dif- 
ferent times. The parameters in Eq. (4.5) for the initial 
state a.t t = 0.0 are Tx = —4.0, ry = 0.0, = 7.0, 
ky = 0.0 and = 0.25. It is seen that the energy 
density profile of the single excitation field which has a 
Gaussian distribution in Fourier space is also Gaussian. 
At a later time t = 6.0 the photon has propagated to the 
right and the intensity profile has spread a little. The 
time evolution is qualitatively the same as obtained in 
the paper S using plane wave quantization. 
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FIG. 10. Time evolution of the energy density distribution 
of a Gaussian one photon state with parameters ~ —4.0, 
Ty = 0.0, k^ = 7.0, ky = 0.0 and A| = 0.25. At t = 0.0 
the intensity profile is Gaussian. Later ai t — 6.0 the photon 
has propagated to the right and the intensity profile is not 
perfectly Gaussian. 



On the contrary to plane wave modes the wavelet 
modes are coupled also for the free field. This makes the 
operation of the free field Hamiltonian to the statevec- 
tor, which is needed in the numerical integration, slower. 

However, most of the coupling constants w^^^ are zero. 
Use of this fact makes the operation much faster. The 
function F^l, (x) in the calculation of the coupling coef- 
ficients can be used to save memory. In general fewer 
wavelet mode functions are needed to represent a local- 
ized field state than plane wave mode functions. 



C. Wavelet mode function initial state 



Next we study the time evolution of the field state 
which at time t = 0.0 has only one wavelet mode excited. 
The mode which is initially excited has parameters s = 1, 
1 = and i = 3. The intensity of the field state is shown 
in Fig. 11. It is clearly localized around r = which is 
expected based on the parameters chosen. The Fourier 
transform for i = 3 mode function is divided into four 



parts at the corners of the frequency interval with a spe- 
cific scaling index s, as is explained earlier and shown in 
Fig. 1^. The intensity profile at time t = 7.5 is shown 
in Fig. 11, The intensity is divided to four main parts 



which all propagate away from the origin. This is under- 
standable based on the Fourier transform of the mode 
function. 
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FIG. 11. Time evolution of a state which at time t — has 
only one Meyer wavelet mode excited. The parameters of the 
mode function are s = 1, Ix = ly = 0.0 and i = 3. The in- 
tensity profile at t = 0.0 is localized at the origin. At t = 7.5 
the intensity profile has propagated mainly to four directions 
away from the center. 



D. Decay of a two level atom 



Next we couple a two level atom to the two dimensional 
field. The interaction Hamiltonian in dipole approxima- 
tion can be written as 



Hi = -T)- E(ro), 



(4.8) 
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where D is the electric dipole moment operator of the 
two level atom 



D 



(4.9) 



and rg the position of the atom. We take the direction 
of the dipole operator to be in the z direction. Using the 
wavelet expansion of the field (3.S) we get 



(4.10) 

In the Hamiltonian the rotating wave approximation 
(RWA) has been used. The approximation can be done 
if the mode functions used are well localized in Fourier 
space. This is the case with Meyer wavelets. The scale 
parameter s determines the frequency interval where the 
Fourier transform of the mode function is nonzero. Typ- 
ically the atom interacts with mode functions of only 
a few, maybe one or two, scales. The mode functions 
are also spatially localized. Because the interaction is 
proportional to the electric mode function evaluated at 
a position of the atom, only mode functions which are 
centered close to the atom interact with it. In this re- 
spect the situation is different compared to the plane 
wave mode functions which are delocalized. In general 
the atom is coupled to fewer mode functions than in the 
plane wave quantization. 

Fig. |l^ shows the logarithm of the excitation prob- 
ability of the atom as function of time. The resonance 
frequency of the atom is w = 10.0 and the dipole cou- 
pling D=0.06. The atom decays energy to the wavelet 
modes. The decay is clearly exponential with a decay 
constant F = 0.18. This corresponds to the theoretical 
value r = \D'^uj'^. 



V. CONCLUSION 



In this paper we have shown how wavelets can be used 
as basis functions in canonical quantization. Different 
mode functions for electric and magnetic fields, which 
are localized both in real and Fourier space, are obtained. 
Mode functions as well as new operators in wavelet basis 
are linear transforms of plane wave mode functions and 
operators. The new annihilation and creation operators 
satisfy bosonic commutation relations. Because the for- 
malism remains the same, it is easy to change the basis 
from plane waves to wavelets. We have applied the the- 
ory to a few example simulations and showed that the 
new basis gives well known results in all cases. In this 
paper we have used wavelet basis for bosonic operators in 
canonical quantization. The same methods can be used 
to change basis also for fermionic operators. A localized 
basis is beneficial for many solid state and semiconductor 
physics problems. 

There are several generalizations and improvements 
of the theory described in this paper. Complex and 
biorthogonal wavelets have some benefits compared 
to real and orthonormal wavelets, for example they can 
be symmetric. One generalization is to use wavelet pack- 
ets or multiwavelets instead of wavelets. Finally it would 
be interesting to compare characteristics of mode func- 
tions of different wavelets. It is also possible to construct 
new wavelets, which have desirable properties, for differ- 
ent problems. 
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